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Migraine is a common disabling headache disorder characterized by recurrent episodes sometimes preceded 
or accompanied by focal neurological symptoms called aura. The relation between two subtypes, migraine 
without aura (MWoA) and migraine with aura (MWA), is explored with the aim to identify targets for 
neuromodulation techniques. To this end, a dynamically regulated control system is schematically reduced 
to a network of the trigeminal nerve, which innervates the cranial circulation, an associated descending 
modulatory network of brainstem nuclei, and parasympathetic vasomotor efferents. This extends the idea of 
a migraine generator region in the brainstem to a larger network and is still simple and explicit enough to 
open up possibilities for mathematical modeling in the future. In this study, it is suggested that the migraine 
generator network (MGN) is driven and may therefore respond differently to different spatio-temporal noxious 
input in the migraine subtypes MWA and MWoA. The noxious input is caused by a cortical perturbation 
of homeostasis, known as spreading depression (SD). The MGN might even trigger SD in the first place by 
a failure in vasomotor control. As a consequence, migraine is considered as an inherently dynamical disease 
to which a linear course from upstream to downstream events would not do justice. Minimally invasive and 
noninvasive neuromodulation techniques are briefly reviewed and their rational is discussed in the context of 
the proposed mechanism. 



Migraine is characterized by recurrent attacks 
of moderate to severe headaches sometimes pre- 
ceded by visual, sensory, motor, or language dis- 
turbances. These so-called migraine aura symp- 
toms are caused by a chemical imbalance in the 
brain that lasts usually not longer than one hour, 
while the longer lasting headaches originate from 
neural activity in the brainstem. Both have 
been observed with non-invasive imaging. In 
this study, it is proposed how these events are 
linked together and that migraine pathophysiol- 
ogy should be considered as a dynamical network 
phenomenon. 



There are several subforms of migraine, yet the two 
main subtypes are migraine without aura (MWoA) and 
migraine with aura (MWA). MWA comprises, for in- 
stance, six subforms, one even without headache. While 
for most migraine sufferers, the head pain is the key man- 
ifestation of their disorder, the migraine aura is also a 
distinctive feature: often visual hallucinations, language 
problems, motor weakness and other short-duration (5- 
60min) changes to sensory modalities or cognitive func- 
tions. These symptoms occur only in MWA, i.e., in about 
30% of the cases and usually precede or sometimes ac- 
company the headache^. Why not in all attacks? This 
raises a related more general question: Is the pain mech- 
anism in MWoA and MWA the same? If so, what drives 
the activation of neuronal pain pathways and how is this 
related to the cause of aura symptoms? And what de- 
termines which sensory modality or cognitive function 
is affected? These interlinked questions about migraine 
pathophysiology remain open or are only unsatisfactorily 
answered in current migraine theories. 

The migraine generator and the spreading depression 
(SD) theory of migraine both focus on distinct patho- 
physiologic events, yet they are not mutually exclusive 
theories of the cause of episodic migraine. The issues 
raised above might be resolved, if both theories can be 
unified. 

From a mathematical point of view, migraine patho- 
physiology involves sudden dynamical transitions, be- 
cause although migraine is a chronic neurological disor- 
der it is characterized by recurrent episodes. In these 
transitions both temporal rhythms and spatio-temporal 
patterns change, which allows only one conclusion: mi- 
graine is a dynamical disease; a term coined to identify 

diseases that occur due to an abrupt change, usually a 
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I. INTRODUCTION 

Migraine is one of the most prevalent neurological dis- 
orders with a life-time prevalence of about 14%. Mi- 
graines can cause substantial levels of disability, ranking 
on a disability scale from 0.0-1.0 at 0.7^. A migraine 
attack can have up to four distinct stages: a prodromal 
phase with difficulty concentrating, yawning, and fatigue; 
an aura phase with all kinds of sensory disturbances; the 
headache phase, often unilateral and throbbing; and fi- 
nally a postdrome of tiredness, difficulty concentrating, 
persistence of sensitivity to light and noise. Yet, despite 
a distinctive clinical picture, migraine continues to be 
underdiagnosed and undertreated. It is estimated that 
the cost for the US and European economies sum up to 
US$19.6 billion and €27 billion a year, respectively. 
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fined as a situation where sudden changes bring the sys- 
tem through a bifurcation 3 . But in a wider sense, if we 
have — or aim for — an understanding of the pathophysi- 
ological mechanism in terms of equations of motion and 
low dimensional phase space structures, we may call it 
a dynamical disease. The sudden change may also origi- 
nate from other dynamical phenomena than bifurcations, 
such as intermittency, noise driven excitation behavior, 
etc.. These various points that lead to sudden transitions 
may also be summarized as "tipping points" . 

As a consequence, migraine pathophysiology needs to 
be considered as a problem at the interface of clinical 
neurology and applied nonlinear science. To this end, 
the migraine generator theory must be extended from a 
theory of brainstem regions as the origin of attacks to 
a network concept where the origin is not localized but 
a dynamical transition into dysfunctional controP^. In 
Sec. |TlJ an extended but still reduced migraine genera- 
tor network (MGN) will be introduced. This network 
is connected to the cortex, where the spatio-temporal 
patterns of SD occur. SD itself can be describes by a 
macroscopic continuous limit of cortical tissue in terms 
of reaction-diffusion systems (Sec. Ill I. How do these two 
dynamical systems interact in migraine, that is, how is 
the central pattern generator network coupled to the cor- 
tical reaction-diffusion system? 

A possible relation is proposed in this theoretical study. 
The pivotal role lies in the driving input of the MGN. 
The MGN is a system that compromises two physiologi- 
cal control subsystems, the trigeminovascular system for 
vasomotor control and associated descending modulatory 
brainstem system for pain control. The vasomotor sys- 
tem may even lead to the ignition of cortical SD in the 
first place, which then in turn drives the MGN. We do 
not consider the cause of the ignition of SD but suggest 
SD — once ignited — drives the MGN and this systems re- 
sponds differently to different spatio-temporal signatures 
of SD, that is, characteristic noxious action of SD leads 
to the subforms MWA and MWoA. In this study, only 
the driving mechanism is considered in some detail in 
Sec. Hn] 

A mathematical model of this dynamical response is 
beyond the scope here, though the MGN itself is pre- 
sented as a reduced scheme explicit enough to support the 
conclusions for neuromodulation techniques (Sec. IV I. In 
general, this dynamical response is a form of central sensi- 
tization, which refers to a transition in a pain network — 
part of the central nervous system — after which it re- 
sponds with increased sensitivity to noxious stimuli (hy- 
peralgesia) and even non-noxious stimuli (allodynia). In 
episodic migraine, the central sensitization is probably 
one of second-order neurons^. Chronic migraine, in con- 
trast, central sensitization is probably a dynamical tran- 
sition in the pain matrix (see below), which also leads 
to an enhanced response to peripheral input but is not 
considered here. 

In Sec. |IV[ minimally invasive and noninvasive neuro- 
modulation techniques are briefly reviewed. These tech- 




FIG. 1. (color online) Networks related to pain, (a) Cra- 
nial circulation (CC): connected by first-order perivascular 
pain-sensitive neurons via the trigenimal ganglion (TG) to 
the trigeminocervical complex (TCC, see b) and receives also 
input from the superior salivatory nucleus (SSN). (b) Brain- 
stem: trigeminal ganglion (TG); trigeminal nucleus caudalis 
(TNC) and C1-C2 regions compromise the trigeminocervi- 
cal complex (TCC); rostral ventromedial medulla (RVM), 
locus coeruleus (LC), ventrolateral periaqueductal grey (vl- 
PAG). (c) Pain matrix^ 2 : thalamus (TH); anterior cingulate 
cortex (ACC); amygdala (Amyg); (PFC); primary (SI) and 
secondary (S2) somatosensory cortices; supplementary motor 
area (SMA); prefrontal cortex (PPC); insula cortex (Insula). 



niques target neural structures in the peripheral and cen- 
tral nervous system that have been introduced in Sec. [TT] 
and Sec. |III| The aim of these techniques is to abort 
migraine. I discuss their rational in the context of the 
proposed driving mechanism of the MGN by SD dynam- 
ics. 

In conclusion, noninvasive techniques for episodic mi- 
graine should be most effective in transcranial stimula- 
tion with three distinct modes of action, for (i) the pro- 
dromal phase, (ii) the ignition and (iii) the acute phase 
of migraine aura. Furthermore, one mode of action in 
transcutaneous stimulation (iv) in the acute pain phase 
of the episodic attacks. Further quantitative approaches 
are needed to find optimized stimulation protocols for 
these situations based on the proposed mechanism. 



II. MIGRAINE GENERATOR NETWORK 

In several respects, pain is different from other sen- 
sations. For instance, normal sensory pathways take a 
chain of three neurons to get from receptors, e.g., sen- 
sitive mechanoreceptors, to the primary sensory cortical 
area. In pain traffic, there are more intermediate nuclei 
and a cascade of descending control mechanisms is essen- 
tial for defining the pain experience^. Therefore a simple 
order of neurons on the way to the cortex is questionable. 
Furthermore, instead of a single, i.e., spatially confined, 
primary cortical area as the first target in the cortex, 
a neural correlate of pain perception seems to be found 
in spatially segregated activity patterns of a widespread 
network of cortical areas. 

Whether there is a network specific for pain percep- 
tion, the "pain matrix" — a term that usually also in- 
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eludes diencephalic structures in the forebrain such as 
the thalamus, see Fig. [I] — or not and instead this net- 
work actually is a largely unspecific sensory "neuroma- 
trix", is still debatecP^. In any case, this matrix seems 
not to be an optimal target in episodic migraine treat- 
ment, e.g., for neuromodulation techniques, because its 
activity patterns are probably truly downstream events 
(but cf. chronification of pain, i.e., >15 headache days 
per month over a 3 month period, at least 8 migrainous, 
and absence of medication overuse). 

The question of the most upstream event is more diffi- 
cult (see Discussion) , and it may only make sense to talk 
about upstream events in the form of transitions in dy- 
namical states of the migraine generator and spreading 
depression (Sec. III). 

While it is worth illustrating all nuclei of this network 
in their respective anatomical place (in approximation), 
it is still necessary to identify a schematically reduced 
network of these nuclei, their neural subpopulations, and 
a vasomotor loop that is explicit enough to convey the 
regulatory control in this network structure amenable for 
mathematical analysis, see Fig. [5] For this reason, the 
following brain regions and their respective functions are 
briefly introduced. 



cal spinal cord, together known as the trigeminocervical 
complex (TCC), see Fig. [I] 

Beyond the TCC, nociceptive traffic ascends by acti- 
vated second-order neurons to further brainstem nuclei: 
the rostral ventromedial medulla (RVM) including the 
nucleus raphe magnus (NRM) in the medulla; the lo- 
cus coeruleus (LC) in the dorsolateral pons (both hind- 
brain); and the ventrolateral periaqueductal grey (vl- 
PAG) in the tegmentum (midbrain). Furthermore, the 
TCC also ascends nociceptive information directly to the 
hypothalamus (HY) and thalamus (TH) with second- 
order neurons along the trigeminohypothalamic tract and 
the trigeminothalamic tract (also know as quintothalamic 
tract), respectively. TH and HY are structures in the di- 
encephalon of the forebrain from where third-order neu- 
rons ascend to the cortex (CTX), the final telencephalic 
station within the forebrain. 

For the sake of completeness, one other pathway of the 
trigeminovascular systems should be mentioned. There 
is a peripheral vasomotor connection from efferents of the 
superior cervical spinal cord to the cranial circulation^!. 



B. Cortical and diencephalic influence 



A. The trigeminovascular system 

The trigeminovascular system, which is key to under- 
standing the network structure of the MGN, is a largely 
bidirectional system. Notwithstanding, it is natural to 
consider as a starting point, if only for convenience, the 
nociceptor and the afferent projection of the first-order 
neuron, that is, the innervation of both dura mater and 
large intracranial vessels by the ophthalmic division of 
the trigeminal nerve. This innervation of the cranial cir- 
culation (CC) terminates in the central trigeminal nu- 
cleus caudalis (TNC) and the C1-C2 regions of the cervi- 
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FIG. 2. (color online) Schematic diagram of pathways in the 
migraine generator network (MGN^ptHIil] 



Before we come to the descending modulatory network 
of the brainstem, let us briefly mention descending path- 
ways from the pain matrix. Cortical projections to the 
periaqueductal grey (PAG) arise from the medial network 
of the prefrontal cortex (PFC) or cortical areas that are 
closely related to this networlP^. The interest here is 
obviously the part which belongs to the MGN, in partic- 
ular the ventrolateral (vl) PAG. The vIPAG evokes pas- 
sive coping reactions (withdrawal, quiescence)^ consis- 
tent with behavior exhibited by migraine patients. Since 
such coping behavior is also readily elicited in animals 
with the forebrain being disconnected from the brain 
sterrP^, a less important role of these descending path- 
ways for the dynamics of the MGN is expected. 

Therefore, we do not include the detailed areas or even 
a closed feedback loop through the pain matrix in the 
attempt to construct a basic structure of the MGN in 
Fig. [2j though this input is indicated by a dotted grey 
arrow. To support this further, it could be shown with 
non-invasive imaging that brainstem activation in hu- 
mans persists in migraine even after pain relief and re- 
duction of symptoms by acute treatmentP. However, it 
might be that not PAG but nearby nuclei in the dorso- 
lateral pons are causing the activation^. 

Notwithstanding, the thalamus is certainly key for 
integration of nociceptive inputs in migraine and the 
pain matrix for pain sensation^. In particular to un- 
derstand chronification of pain, closed feedback loops 
through these areas are probably important and these 
diencephalic or cortical areas may even primarily be the 
"upstream" target areas for a therapeutic rational in 
chronic migraine. 
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C. Descending modulatory network 

Given that PAG is an area to which various cortical 
and diencephalic structures project, it is important to 
note that PAG does not directly project itself to the 
TCC. The descending projection from PAG modulate 
ON and OFF neurons in RVM. As indicated by the name, 
these ON and OFF neurons in the PAG-RVM descending 
control system can be facilitatory and inhibitory depend- 
ing on how the balance shifts in either direction^. A 
descending facilitation of TCC can provide a mechanism 
for central sensitisation. The development of cutaneous 
allodynia suggests that the balance shifts in favor of fa- 
cilitation explaining the increased sensitivity of the skin 
to non-noxious stimuli in migraine^. 

It has been questioned whether the PAG or rather the 
nuclei in the dorsolateral pons, like LC, are responsible 
in this descending modulatory network and it was sug- 
gested that "[t]o advance the brainstem migraine gen- 
erator theory from the opinion phase to being evidence 
based, answers should be provided to questions as: (a) 
How would increased activation in the PAG drive or pro- 
duce migraine?"^. Quantitative approaches in a form of 
a mathematical model similar to the gate theory of first- 
order neurons^ may contribute to this question. The 
important notion stressed here is that the extended mi- 
graine generator network theory is essentially a dynamic 
network theory. The MGN is a central pattern gener- 
ator, that is, a neural network that produces rhythmic 
patterned output perceived as pain without sensory feed- 
back. However, such feedback is still needed, as proposed 
here, to drive the transition into central sensitization. 
The explanatory power of a mathematical approach lies 
in the predictive power, but it certainly also needs the 
clinical testing to be evidence based. 

D. Vasomotor control 

The key vasomotor control is probably played by a 
parasympathetic cluster of neuronal cell bodies located 
in a fossa (ditch) in the skull, the sphenopalatine 
ganglion (SPG). The simplest putative structure of the 
MGN (Fig. [2} is completed by the SPG and its second- 
order preganglionic neurons in the superior salivatory 
nucleus (SSN). It was suggested that migraine triggers 
typical for the prodromal phase either activate or even 
originate in a number of brain areas whose projections 
converge on the SSN (dotted grey arrow) and which are 
functionally positioned to produce migraine symptoms 8 . 



III. SPREADING DEPRESSION THEORY OF 
MIGRAINE 

SD is essentially a slow (about 3mm/min = 50/im/sec) 
reaction-diffusion wave in gray matter tissue. In the cor- 



tex, SD is accompanied by a pronounced hemodynamic 
response of increased regional cerebral blood flow (hyper- 
mia) for about 2min and a long lasting, ^2h, decrease 
(oligemia^. 

The cortical tissue SD traverses is massively perturbed 
in its ion homeostasis 19 . If the ion flow across the mem- 
brane in this region is taken to estimate relative changes 
in Gibbs free energy of the tissue, SD reveals itself as 
"a twilight state close to death"^, i. e., during the peak 
of SD the cortical state is similar to the one in stroke 
or ischemia. The ion flow causing this state is far more 
dramatic than in any other neurological disorders, for in- 
stance, compared to the ion flow during ictal epileptic 
activity. 

However, the cortex is not pain sensitive. It was 
therefore sugg ested that SD may trigger the pain phase 
indirectly^al. 

A. SD and pain pathways 

In several animal models, possible pathways were in- 
vestigated to see how SD might cause pain by activating 
nociceptors in sensitive intracranial tissues and subse- 
quently activate the TCO^S. In short, it was suggested 
that during SD the blood-brain barrier becomes more 
permeable, allowing glutamate, potassium ions, hydro- 
gen ions, nitric oxide and other noxious substances or in- 
flammatory mediators to diffuse from the surface of the 
cortex into the meninges where they activate nocicep- 
tors of the cranial circulation. Although this mechanism 
has been criticized^, the animal models of migraine show 
that SD activat es the first- order neurons in the trigeminal 
nucleus caudali^- 21 * 22 * 25 * 26 ! In a newly published article, 
it was found that SD can trigger headache by activating 
neuronal Panxl channels, a megachannel that is a struc- 
tural component of gap junctions^. It is furthermore 
discussed, whether nuclei in the MGN in turn can cause 
spreading depression in the cortex. For a summary, see 
RefP. 



B. SD, aura, and silent courses 

Even the apparently solved question of how SD relates 
to the aura phase has still some merit. It is clear that the 
aura is caused by SD. At least after successfully imaging 
the blood flow in migraine aura, it became indisputable 
that SD is the electrophysiological correlate of migraine 
aur£ (29H3i] i n fact, already in 1945, Leao wrote: "The lat- 
ter disease [migraine] with the marked dilatation of major 
blood vessels and the slow march of scotomata [sensory 
blindness] in the visual or somatic sensory sphere is sug- 
gestively similar to the experimental phenomenon here 
described [SD] , in spite of the fact that known scotomata 
are still felt to be vasoconstrictor in nature"^. It was 
not known in 1945 that in migraine, a vasoconstrictor 
response (oligemia) is also present. 
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FIG. 3. (color online) Paradigms of SD pattern in migraine. 
Narrow band of excitation (red), refractory zone (yellow), 
disperse neurovascular feedback (green), (a) Typical text- 
book illustration: SD engulfing large parts of one hemisphere; 
adopt ed fro m RefP. (b) Localized SD wave in a full-scale 
attadPEH. 



To date, what remains disputed is whether SD is even 
present but clinically silent in the migraine subtype 
MWoA. This is referred to by the contested notion of 
"silent aura'^. If true, it is SD which stays silent, the 
aura is simply not present. Therefore "silent SD" is a 
term less perplexing. Anyway, such silent SD courses are 
supported by (a) a well-documented case of blood-flow 
changes that were likely the result of SD observed in a 
spontaneous migraine headache without auraP^ and (b) 
the fact that SD can cause pain and MWA and MWoA 
are likely to share the same pain mechanism. 



To resolve these questions about the relation of SD 
to the pain and aura phase, we recently suggested that 
which of the two migraine subtypes develops, i. e., MWoA 
or MWA, is dependent upon the spatio-temporal pattern 
(shape, size, and duration) of SEW This does not rule 
out the possibility that SD itself is triggered by events 
in the brainstem and therefore we contribute only little 
to the questions of the most upstream event. However, 
the crucial part missed in the current controversy might 
be the spatio-temporal pattern of the cortical tissue be- 
ing recruited into a depleted state. All too often, it is 
completely ignored that SD usually propagates in the 
human gyrified cortex as a discontinuous wave segment 
as observed with fMRI in migraine and reported by pa- 
tients as visual field defect pSEUSID see pjg [JJ-, j n am _ 
mal models, SD is engulfing the whole lissencephalic cere- 
bral hemisphere but already in gyrencephalic cat brain 
only the first SD waves do this and succeeding secondary 
SD waves, which propagate in relative refractory, i.e., 
less susceptible, cortex, often remained within a confined 
origirPZl. 



C. Localized SD waves 

Both fMRI data and reports on visual field 
defects^TMM] clearly show that SD is not only spatially 
confined in humans, but SD is a discontinuous wave — a 
particle-like wave. A particle-like wave is a wave front 
with two open ends, sometimes also called dissipative 
solitorP^l. The existence of such localized waves is a 
change in the pattern forming paradigm indicating a fast 
diffusing second inhibitor related to the Turing bifurca- 
tion. 

From a clinical perspective, migraine aura symptoms 
are focal symptoms. If SD were to spread out in all di- 
rections, the affected cortical area would increase quickly 
and symptoms should become more generalized. Local- 
ization does not merely restrict the symptoms caused by 
SD; but localization also allows SD to assume a much 
larger variety of spatio-temporal patterns compared to a 
continuous circular wave front. This is reflected in the va- 
riety of the courses of focal aura symptoms^ and maybe 
also in the fact that SD could stay silent under certain 
conditions on the localised pattern (Sec. Ill B,F) and, 
furthermore, that SD does not cause pain under other 
conditions (Sec. Ill F-G). 

For these reasons, it is an essential feature that SD 
stays localised even when it propagates a long distance 
in one direction. In contrast, if SD spread out in all 
directions, as to date shown in all modern migraine text- 
books (published or revised in the last 5 years) , the shape 
would essentially be circular (cortical folding has some in- 
fluence) and, furthermore, size and duration both would 
depend on a single parameter, the radial distance SD 
travels. The course that is usually sketched starts from 
the occipital pole of the cortex and SD propagates in 
a full-scale MWA attack up until the frontal lobe, see 
Fig. |3^. If SD invaded the complete occipital, temporal, 
and parietal lobe, this would not match with the reported 
courses of neurological symptoms during the aura phase, 
not even in a full-scale attack, cf. RefP^ and it is also 
implausible that only a subset of the nearly completely 
depleted cortical tissue results in some form of clinical 
manifestation. 

A simple calculation illustrates the different scaling be- 
havior of events. In a full-scale MWA attack, the aura 
phase can last up to one hour. This corresponds to 
a length I of li = 18cm cortex being paced. Even if 
we halve this distance I2 — 9cm — either because SD 
is slower or lasting only 30min and symptoms reverber- 
ate for a full 60min — this would still correspond to 20% 
(80% for li = 18cm) of the surface area of one cerebral 
hemisphere. In contrast, only 1.5% (3% for l\ — 18cm) 
of the surface are affected in a full-scale MWA, if the 
particle-like SD wave is on average 2cm wide. Note that 
the progressive development of aura symptoms is, except 
for the ignition phase of SD, probably mostly due to the 
fact that SD starts in regions of high cortical sensory 
magnification and travels into regions with lower cortical 
magnification^]. 
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D. Generic reaction-diffusion model 

In accordance with the noninvasively imaged SD pro- 
gression and reported visual field defects, we proposed 
a model in which in each attack a particular pattern is 
formed by a discontinuous wave segment that spreads 
out only in one direction determined by the initial 
conditions^ and guided by the folding pattern (unpub- 
lished results). The ignition patterns that trigger these 
simulated attacks are given by an assumed hyperactiv- 
ity in a cortical feature map of the visual cortex used as 
initial conditions for the model. These feature map pat- 
terns were also suggested to be hyperactive during SD 
causing visual hallucinatory percepts 4 ^. 

In this generic reaction-diffusion mechanism of ex- 
citable media, the resulting pattern forming processes 
following these initial conditions are described in ab- 
stract terms of activator-inhibitor kinetics of two vari- 
ables u and v, respectively. In fact, simple activator ki- 
netics of u date back to a mechanism of SD described in 
1963 by Grafstein, Hodgkin and Huxley (GHH)P. The 
GHH model assumed u = [K + ] e , i.e., the extracellular 
potassium ion concentration to be the activator u. We 
extended this scheme (for details of the full model see 
RefP): 

— =u- -u 3 -v + DV 2 u, 1 
at 3 

~ =e (u + P + kJh(u) dxdy) , (2) 

with D the diffusion coefficient, two parameters e (time 
scale separation) and (3 (threshold), H being the Heav- 
iside step function (discontinuous function whose value 
for negative (positive) arguments is (1) and 0.5 for the 
argument zero). The integral term is therefore a surface 
area measure of the medium being in a state with u larger 
than zero that can globally increase the threshold. 

Both activator u and inhibitor v are lump variables. 
Unlike in the original GHH model, it is not necessary or 
even possible to identify particular physiological quan- 
tities, like [if + ] e , with them. Instead, the activator u 
should be viewed as the bistable energy state. One stable 
fixed point corresponds to the maintenance of homeosta- 
sis far from thermodynamic equilibrium and the other 
stable fixed point to the state where the cellular Nernst 
reversal potential are depleted (thermodynamic equilib- 
rium). In the presence of an inhibitor v, which was not 
included in the GHH model, the second stable fixed point 
becomes a transient state. In other words, the inhibitor 
v is a recovery variable. It is related to ion pumps that 
drive the tissue far from thermodynamic equilibrium and 
(re)charge Nernst reversal potentials. We set the parame- 
ters as follows: time scale separation e — 0.04, the thresh- 
old /3 = 1.32, and the mean field coupling K — 0.0 03. 

The details of the model are discussed elsewher e! 34 * 35 !. 
In short, for K = 0, the bistable activator kinetics 



with diffusion and the linear inhibitor kinetics (immo- 
bilised) provide the simplest reaction-diffusion model of 
activator-inhibitor type for the cortex as an excitable 
media exhibiting patterns as seen during SD in animal 
models. The homogeneous steady state is a stable so- 
lution representing the healthy cortical state. A critical 
neuronal mass is needed to start ignition and support 
sustained propagation of a Gibbs free energy-depleted 
state, which is associated with the state of SEBSl. F or 
K = 0, this depleted SD state would engulf all the tissue 
after an initial ignition. The suggested pattern forming 
mechanism of localized SD is closely related to the criti- 
cal mass needed for ignition, called a critical nucleation 
solution, see Fig. [4^,. 

This critical nucleation is controlled for if ^ by a 
widely spread out inhibitory feedback, Fig. [4]d,c. This 
inhibitory feedback is approximated by a mean field term 
K J H (u) dxdy, a term proportional to the depleted sur- 
face area region in the cortex. The physiological signifi- 
cance of this last term in the rate function of the lump 
variable v is associated with the neurovascular feedback 
(e.g ., as observed during hypermia of increased blood 
fk W 29 * 31 *) in SD. For the sake of simplicity, we chose it to 
be proportional to the surface area of the depleted state, 
although the neurovascular feedback is mainly driven by 
the neuronal hyperactivity in the rise of the SD front. 
Therefore, a measure that is proportional to the length 
of the SD front may be more appropriate. However, we 
tested K J H (u) H (1 — v) dxdy, which is an approxi- 
mated measure of the length of the SD front — the ac- 
tivator must be in the excited state while the inhibitor 
is still low, as characteristic for the hyperactivity in the 




FIG. 4. (color online) Schematic orbit structures in phase 
spaced see text, (a) Caricature of excitable media, (b) Basin 
of attraction change when local excitabilty parameter change, 
(c) Fast secondary inhibition creates localized structures, (d) 
Ghost behaviour when traveling wave solution (node) collides 
with nucleation (saddle): the two trajectories are associated 
with migraine with aura (MWA) and with migraine with typ- 
ical aura without headache (MxWA). (e) Direction of cortical 
area affected by SD. (f) Trajectory of a circular wave associ- 
ated with migraine without aura (MWoA). 
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rising front — with similar results. 



E. Transient waves and slow dynamics 

Changing a parameter of this model without the mean 
field inhibitory feedback (K = 0), essentially varies the 
basins of attraction of the existing solutions, e.g., mak- 
ing the medium less susceptible by increasing the basin 
of attraction of the homogeneous steady state, Fig. [4]b. 
For K ^ 0, the traveling wave solution itself is changed 
such that the traveling wave solution becomes localized, 
as described in the previous section, Fig. |4j:. For suf- 
ficiently large K, this stable localized solution collides 
with its nucleation solution (saddle), which leads to well 
known phenomenon of slow dynamics named ghost be- 
havior after a saddle-node bifurcation, Fig. |3Ji. This is 
the origin of the transient nature of the SD signatures in 
Fig. [5] that we propose to lead to migraine without aura, 
migraine with aura, and migraine with typical aura and 
without headache (MWoA, MWA, and MxWA, respec- 
tively) . 

In particular, MWA, and MxWA correspond to trajec- 
tories that pass the ghost of the slow manifold (lightblue 
in Fig. [4jd) . In contrast, solutions that are shorter lasting 
but with larger MIA (cf. Fig. J4J?) do not pass the ghost 
of the slow manifold Fig. [4]:; they correspond to MWoA. 

In Sec. IV D, this is further discussed in the context of 
the therapeutic strategy for neuromodulation. 

The macroscopic features of SD are well described by 
this model. The affected cortical area will not grow in 
time and SD stays localized while propagating (for fur- 
ther discussion of the usefulness of such a toy model, see 
Sec. IV D |. It is worth noting that it is not without irony 



that both terms in the name of SD are rather misleading: 
the important neural activity is a hyperactivity that only 
in the depleted phase becomes "depressed" , and this cor- 
tical state is only in the ignition phase "spreading" , the 
sustained main process is a traveling wave. 



F. Link to migraine subforms 

We performed statistical analysis of the spatio- 
temporal development governed by the generic reaction- 
diffusion model in Eq. (l)-(2) that followed 8000 different 
natural initial perturbations of local hyperactivity as SD 
triggers. These perturbations were introduced by an in- 
crease of the resting value of u. Based on that, we predict 
that certain features related to shape, size, and duration 
of SD (Fig. [5^) determine the aura phase while others 
determine the pain phase in migraine. 

It is important for further discussion that we firstly 
divide the subtype MWA into its two subforms "typi- 
cal aura with migraine headache (1.2.1, in the interna- 
tional classification)" and "typical aura without headache 
(1.2.3)", which we refer to as MxWA. Therefore, in the 
remainder of the text, we refer to 1.2.1 as MWA, i.e., 
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FIG. 5. (color online) (a) Spatio-temporal signatures of SD as 
obtained from a reaction-diffusion model with mean field cou- 
pling, Eqs. (l)-(2); cf. Fig. 4E in RefP. The model param- 
eters are chosen such that only transient waves exisfl^, note 
that we have scaled the dimensionless model such that the 
units of time are roughly corresponding to minutes, and sur- 
face area to cm 2 and the wave spreads with about 3mm/min. 
Four stereotypical courses are depicted following a local per- 
turbation of a homogeneous steady state: (1) subthreshold, 
SD dies out quickly without initial spread; (2) superthresh- 
old perturbation, SD dies out after a few minutes, if its front 
does not break open, corresponding to migraine without aura 
(MWoA); (3)-(4) superthreshold perturbation, SD propagates 
for up to 30min (~15 in (3), ~25min in (4)), if the SD front 
breaks open, corresponding to migraine with aura (MWA). 
(b) Statistical analysis of 8000 events plotted over MIA and 
TAA, see text. 



we explicitly assume the presence of headache. MxWA 
is estimated to account for about 5% and MWA for 25% 
of the 30% of the subtype migraine with aura (1.2). The 
number of unreported cases of MxWA might be quite 
large, though. 

Together with MWoA, three combinations exist. We 
predict that each migraine subform has a specific sig- 
nature of SD. Characteristic for this signature is not 
only the spatio-temporal pattern of the course of SD in 
the gyrified human cortex but also the hemodynamic re- 
sponse and the amount of noxious substances that diffuse 
into the meninges. The cortical folding pattern should 
also influence this signature. We are currently perform- 
ing a detailed analysis of this influence. 

In particular, we predict that a sufficiently large sur- 
face area must be instantaneously affected and thus de- 
pleted by SD to lead into the pain phase in migraine. 
Therefore, we propose that if the maximal instanta- 
neously affected area (MIA, Fig. [5]d) is too small, the 
cascade of subsequent events causing sustained activa- 
tion of trigeminal afferents is not initiated. The ratio- 
nal behind this is that the transmission of substances 
(noxious or inflammatory mediators) in the direction per- 
pendicular to the cortical surface into the pain sensitive 
meninges should be significantly convergent to reach nox- 
ious threshold concentration and initiate central sensiti- 
sation of second order neurons, Fig. [6] While the flow 
driven by small MIA is sufficiently diluted and therefore 
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FIG. 6. (color online) Schematic representation of a corti- 
cal cross section with meninges and skull. Top: large MIA; 
bottom: lower MIA, not to scale; cortical thickness is at the 
occipital lobe about 2mm; diameter of MIA pain threshold is 
approximately 1.78cm. 



tolerated. This makes a prediction that can be tested by 
noninvasive imaging. 

Furthermore, only if SD assumes a spatio-temporal 
form that is lasting long enough (>5min) and therefore 
also propagating farther (>1.5cm), SD will cause notice- 
ably aura symptoms. In fact, this is not a prediction. 
It merely reflects the diagnostic criteria of migraine aura 
given by the International Classification of Headache Dis- 
orders: "focal neurological symptoms that usually de- 
velop gradually over 5-20 minutes" . So any neurological 
events that last less than 5min are usually not diagnosed 
as migraine aura. 

The pain phase is completely suppressed in cases of 
MxWA, if the ignition is not causing a large MIA. 
Such a temporal evolution in Eq. (l)-(2) is expected for 
initial perturbations that are just large enough to be 
suprathrcshold stimulations and at the same time are 
asymmetric to lead directly into a discontinuous wave 
(Fig. [4jZ> ) that is slowed down. As a consequence from 
the model, MxWA are particularly long lasting auras, 
again a prediction that can be tested. 



G. Pain threshold determined by SD reaction-diffusion 
model 

The distribution of spatio-temporal SD patterns that 
results from local perturbation led us to suggest a pain 
threshold determined by MIA. This value can be esti- 
mated to be about MIAss 2.5cm 2 . In other words, if SD 
is never depleting the transmembrane gradients in cor- 



tical tissue over a surface area of more than 2.5cm 2 (at 
any given time), the pain cascade is not initiated. This 
corresponds to a circular area with a diameter of about 
d ~ 1.78cm or 0.1% of the total cortical surface. The 
course (2) shown in Fig. [5] has a MIA of about 4.6cm 2 , 
while the course (3) and (4) are just above the pain 
threshold. 



In a nondimensionalized generic model, we can only 
give an estimate of the surface area of the pain threshold. 
In fact, we rely on some self-consistency which can have 
biased the result as discussed briefly in the following. 



To estimate surface area in units of cm 2 , the nondimen- 
sionalization of the SD model needs to be considered. We 
could have started with bistable kinetics for the activa- 
tor, similar to that as suggested in the GHH model for 
u = [K + ] e with three roots at 3mM and 55mM, for the 
stable fixed points, resting state and depleted state, re- 
spectively, and lOmM for the threshold, a value that is 
under physiological conditions well pr eserved and there- 
fore called the ceiling level of [if + ]J=*l Also the diffusion 
coefficient for potassium ions in the brain can be used. 
This would result in an equation structurally similar to 
Eq. (1) but in a dimensionalized form. If we were to add 
linear kinetics for an immobilized inhibitor (for the sake 
of argument, set in Eq. (2) K — 0), there is a total of 7 
parameters (4 for the cubic rate function of u, the diffu- 
sion coefficient, and 2 for the linear rate function of v). 
We can scale 4 quantities : the 2 concentrations of u and 
v as well as time and space. Therefore, only 3 free pa- 
rameters remain in a nondimensionalized model. Since 
the rate function in Eq. (2) (still with K — 0) is only a 
function of u and not also v, this reduces further to the 
two parameter /3 and e — for the weak limit, we consider 
here, this is canonical. Weak limit refers to the fact that 
the threshold of the local dynamics (D = 0) is large and 
the bifurcation that determines the oscillatory behavior 
(type I or type II) is irrelevant. 



For any set (s,/3) in the parameter plane, we can re- 
construct dimensionalized space and time scales, because 
we know both the speed of SD and the width of the wave 
profile^. This procedure can be performed in a similar 
way with K ^ 0, although we select by the choice of the 
mean field feedback the size of the critical surface area. 
This is because the choice of K selects a specific path 
in parameter space through the saddle-node bifurcation 
and therefore also the "ghost" behavior, which is indi- 
rectly linked to the obtained value of the pain threshold, 
see RefJS for a detailed discussion. For the values used 
here, we have approximated the dimensionalized time 
and space scales td and Xd (and yd) by td = t/10 and 
Xd = x/8, where t and x are nondimensionalized time 
and space scales from Eqs. (l)-(2), respectively. 
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IV. DYNAMICAL DISEASE AND 
NEUROMODULATION 

It is often not helpful or even meaningful to talk about 
an upstream cause for a composite dynamical system be- 
ing out of control, if the control is established in closed 
loop feedback between the composites. It is likely that 
migraine episodes are caused by events that should be de- 
scribed as transitions in the dynamical state of the brain. 
And it can take various triggers to make this transition 
but the cause is inherently in the dynamics. Therefore, 
let us address the question: Why do we want to under- 
stand migraine as a dynamical disease? 

"The significance of identifying a dynamical disease is 
that it should be possible to develop therapeutic strate- 
gies based on our understanding of dynamics combined 
with manipulations of the physiological parameters back 
into the normal ranges'^. This quote provides a general 
answer to this question and to why therapeutic strate- 
gies that rely on the dynamics will be considered in this 
section. 

The focus in this study is on neuromodulation tech- 
niques. One reason is because pharmacological treat- 
ment is not suitable for any direct dynamical treatment 
at high frequencies because most fast time scales are 
damped due to pharmacokinetics, that is, absorption, 
distribution, and metabolism, before pharmacodynam- 
ics sets in — which of course can act indirectly on fast 
time scales. These issues are bypassed in neuromodula- 
tion. Neuromodulation is also usually (though not al- 
ways, see below) spatially confined. Since dynamical dis- 
eases are caused by transitions in both temporal rhythms 
and spatio-temporal patterns, neuromodulation is the 
natural ansatz for their treatment. "The headache future 
is bright for neuromodulation techniques, if we manage 
to understand how they work." A quote from a slide by 
Jan Schoenen (with permission). 

Another reason is that various minimally invasive and 
noninvasive neuromodulation techniques are available 
and have been tested clinically in migraine^. To date, 
these treatment options are usually only considered in 
migraine in chronic cases of several disabling attacks per 
month that are refractory to other current treatment. 
But they might become available in episodic migraine, if 
safety issues are appropriately resolved. 



A. Phase-dependent mode of action 

Possible target regions are the brain structures de- 
picted in Fig. [2| At least two distinct dynamical situ- 
ations during a migraine episode should be distinguished 
because the optimal stimulation protocol is likely to de- 
pend not only on the mode of action but also on these dif- 
ferent phases. The ignition phase (IP, Fig. [8^) where SD 
is actively controlled by hyperemia, and the longer last- 
ing acute phase (AP, Fig. [8Jd), including oligemia (~ 2h) 
and sustained pain activity, which can last in migraine 



from 4 to 72h. 

The IP lasts about as long as SD circumscribes a large 
surface area in the cortical tissue, which is assumed to 
be above the pain threshold, Fig. [7| During the acute 
aura phase, there is a cortical spatio-temporal pattern 
of both hyperemia and oligemia related to the SD pat- 
tern. This pattern might project even to the nuclei of 
the MGN, that is, it is mapped in a topographic repre- 
sentation. The acute pain phase can overlap and last up 
to 72h. Therefore Fig. [5] oversimplifies these phases but 
gives a first hint on the dominant dynamics. In partic- 
ular, Fig. [8^1 shows the traffic that eventually results in 
central sensitization, while Fig. [SJ^ shows the facilitated 
pain traffic during central sensitization (cf. Sec. EB. As 
a consequence, different neuromodulation targets and/or 
stimulation protocols should be developed and applied 
for IP and AP. Various techniques are available, as briefly 
reviewed in the following. 



B. Invasive neuromodulation 

It is worthwhile to note that patients with chronic pain 
syndromes — but headache-free — developed headaches 
with migrainous features following implantation of deep 
brain stimulation (DB S) ele ctrodes that target in peri- 
aqueductal grey (PAG^^. This was taken as evidence 
for the significant role of PAG in the migraine generator 
(see Sec. [nl cf. Ref.^). Although posterior hypothalamic 
DBS is applied to cluster headache, DBS as an invasive 
method does not play a significant role in headache treat- 
ment to date. The 3% risk of bleeding is only within the 
upper range reported for DBS in movem ent d isorders and 
other than in these disorders (cf. RefsP^Hl), the mode 



ignition phase AP: acute phase 



IP Stim. 
AP Stim. 




FIG. 7. (color online) Simulated course of cortical surface 
area depleted by during SD. Insets: snapshots of simulated 
SD wave in 2D, top view, (a) ignition shortly before MIA 
is reached; (b) ring wave before breaking into discontineous 
wave; (c) discontineous wave; (d) propagation near "ghost" , 
note the new reference frame; (e) propagation behind the 
ghost, new reference frame; (f) collapse. 
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FIG. 8. (color online) Suggested facilitated flow of pain 
traffic in the reduced migraine generator network, (a) Flow 
during ignition of SD; (b) after central sensitisation. 



C. Noninvasive neuromodulation 



Transcutaneous and transcranial stimulation are the 
two noninvasive modes. Transcutaneous electrical nerve 
stimulation (TENS) modulates neural activity by tar- 
geting noninvasively peripheral nerves (cf. ONS). In mi- 
graine, a supraorbital transcutaneous stimulation (STS) 
was tested with a portable device that looks like a sil- 
ver headband and targets the ophthalmic division of the 
trigeminal nerve — the first-order neuron that ends in the 
TCC. The stimulation was tested with biphasic rectan- 
gular AC impulses at 60Hz. The therapeutic gain (26%) 
is within the range of those reported for other preven- 
tive drug and nondrug antimigraine treatments^!. But 
like with ONS, the mode of action or even quantitative 
methods to test paradigms are currently missing, though 
it seems suggestive that if central sensitization of second- 
order neurons are causing the pain, this method is sup- 
pressing not the sensitivity to noxious stimuli or non- 
noxious stimuli but the inbound sensory traffic. 



of action remains unclear for headache. In fact, even 
placebo has not been conclusively ruled ouP^. 

Therefore, current development efforts are primarily 
on minimally invasive, and noninvasive (see below) neu- 
romodulation. In IP, the sphenopalatine ganglion (SPG, 
Sec. II D) plays a central role, therefore minimally inva- 
sive sphenopalatine ganglion stimulation might be suited 
in particular in this phase. An implantable battery-free 
SPG stimulator has been developed to apply on-demand 
stimulation for the treatment of migraine and other pri- 
mary headache^. In a pilot study (n=ll), the potential 
of electrical stimulation of the SPG has been shown, 2 
patients were pain-free within 3 minutes of stimulation, 
3 had pain reduction, 5 had no pain relief, 1 was not 
stimulated^!. However, in this pilot study, the headache 
was allowed to intensify up to 6 hours before stimula- 
tion was initiated, which means that mode of action is to 
be considered in the AP paradigm, where SPG probably 
takes a less direct influence, see Fig. [HJd. 

A promising target is the TCC in both AP and IP. 
With occipital nerve stimulation (ONS), the superior cer- 
vical spinal cord of the TCC is accessible with minimally 
invasive peripheral neuromodulation. Subcutaneously 
(beneath skin) placed electrodes that start at the level of 
CI and t ransv erse to C3 affect the cor respo nding occip- 
ital nervtPIH. ONS is currently tested^ 21 . To develop 
therapeutic strategies depending on a specific mode of 
action in ONS, quantitative approaches for the migraine 
generator have to be further developed to understand 
how the facilitated pain traffic during central sensitiza- 
tion (Fig. [8]d) can be influenced via the ophthalmic divi- 
sion of the trigeminal nerve. 



Transcranial stimulation targets the cortex. Two non- 
invasive neuromodulation techniques are available in mi- 
graine, transcranial magnetic stimulation (TMS) and 
transcranial electrical stimulation (TES). A carefully 
sham-controlled larger trial with a portable TMS devise 
was undertaken as a promising noninvasive neurostim- 
ulator for disrupting migraine attacks^. The simplest 
technology with the least safety issue is probably TES, a 
method that therefore might pave the road for noninva- 
sive portable devices for episodic migraine in less disabled 
patients in the future. Two TES versions exist, transcra- 
nial direct current stimulation (tDCS) and alternating 
current stimulation (tACS). The lack of spatial specificity 
is, however, a critical point. Cathodal tDCS inhibits, 
anodal tDCS increases neuronal firing. Preliminary evi- 
dence was found for patients with chronic migraine hav- 
ing a positive, but delayed, response to tDCS applied to 
motor (anodal) and oribitofrontal (cathode) cortices^!. 
Computational simulations of current flow through brain 
regions were used to interpret the effects. With the help 
of such simulations, a principal understanding of pain 
modulation can be gained, whether this modulation is 
due to the effect of TES on SD or on the dynamics of 
the pain matrix (in chronic pain). Furthermore, the ef- 
fect of TES on deeper structures beyond the immediate 
cortical target regions can be estimated, including the 
cingulate, insula, thalamus, and brainstem--". This is a 
clinical research model that could serve as a guide to 
future investigations. Effective prophylactic therapy in 
migraine was also studied with cathodal tDCS over the 
visual cortex with the anode overlying the motor and 
sensory cortices^. When alternating current (AC) is ap- 
plied over the scalp it is known as tACS, but tACS has 
not been clinically tested in migarine. 
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D. Model-based control 

The above mentioned techniques need a model of the 
pathophysiological events in the target region to optimize 
the stimulation protocol. This is even mandatory in a 
closed loop feedback control, such as a Kalman filter, 
but it is also important for open loop control, which is 
the current mode in all neuromodulation for migraine. 
Simply said, given that we know how the hardware can 
interface with the wetware, what is the best software? 

Quantitative models for SD are already developed in 
great detail. Not only macroscopic models are available, 
such as suggested and briefly introduced here, but also 
detailed models including the electrophsiological events 
on a cellular level, for a review see RefP^l This opens up 
a model-based approach to target SD in the cortex. 

Macroscopic equations of motion for SD, as Eqs. (1)- 
(2) — which also are considered as "toy models" — have 
two advantages. First, they allow extensive computa- 
tional simulations of spatio-temporal pattern that extend 
in the case of SD in migraine over several centimeters and 
last up to one hour. A microscopic model of SD describ- 
ing the dynamics on the millisecond scale of single cells 
to calculate the noxious signatures that are transmitted 
into the meninges during one hour over several square 
centimeters is to crack a nut with a sledgehammer. How- 
ever, the main advantage of such generic models lies in 
the fact that they allow insight in the phase space struc- 
ture of the whole class of models they represent, Fig. [4] 
In this sense they are generic and well suited in particular 
to develop open loop control stimulation protocols. 

Given the predicted spatio-temporal signatures follow- 
ing from the model in Eq. (l)-(2), stimulation protocols 
in the ignition phase and the acute phase have different 
control aims. The acute phase can be further subdivided 
into one with the aura present and one with only the 
headache. Furthermore, the prodromal phase exists. 

The most obvious stimulation protocol can be derived 
for the acute phase with the aura present. In this phase 
the slow dynamics are governed by ghost behavior after 
a saddle-node bifurcation, Fig. [4jl. Noise can speed up 
the transition time in this case and shorten the aura. 
Although this is not the major contribution to the pain 
phase as SD is below the pain threshold, once this thresh- 
old is crossed during ignition, it is likely beneficial to 
suppress any further release of noxious or inflammatory 
mediators by the decaying SD wave. 

According to our prediction, the most critical phase is 
the ignition phase. In this phase the pain threshold will 
be usually crossed. It is also the shortest lasting phase of 
only a couple of minutes. This phase is therefore likely to 
be missed, if neuromodulation is not used in a closed-loop 
control with some sort of feedback. 

Therefore, the prodromal phase becomes in particu- 
lar valuable for prophylactic protocols. In this phase, 
cortical excitability might be temporally decreased to 
lower the upcoming ignition. Prophylactic protocols with 
cathodal transcranial direct current stimulation of the 



visual cortex are currently testecP^l. Similar approaches 
should be tested with TMS, which can be more spatially 
selective to interact with cortical hot spots as targets. 

In the acute pain phase, SD is not a target anymore. 
Therefore the model contributes nothing to this phase, 
except, as already noted above, that according to the 
proposed mechanism transcutaneous stimulation will in 
this case target the inbound sensory traffic. 



V. SUMMARY 

The migraine generator (Sec.|n]) and the spreading de- 
pression (Sec. Ill) theory of migraine are not mutually 
exclusive theories of the cause of episodic migraine. Suf- 
ficient evidence speaks to the fact that both theories can 
be validly maintained. It may be that the migraine sub- 
types MWoA and MWA have different pathophysiological 
mechanisms, though this seems rather unlikely. Likewise, 
it cannot be ruled out that the phenotype of each classi- 
fication subtype has different pathophysiological mecha- 
nisms. 

The proposed mechanism still allows space for either 
of these possibilities. The unified framework, however, 
does not rely on either of these options. We propose that 
stereotype spatio-temporal signatures of SD (Fig. [5| can 
selectively drive rhythms in the MGN, which in turn is 
part of a larger network (Fig. [5]) that controls the spatio- 
temporal dynamics of SD in the first place. 

The question of the most upstream events in episodic 
migraine lead to a long standing controversy. To my 
mind, the settlement of this controversy points strongly 
to the need to complement clinical and experimental 
data with quantitative approaches, in particular com- 
putational simulations and mathematical analysis of the 
dynamical interaction between the migraine generator 
and spreading depression. This study can only outline 
a framework for such quantitative approaches. 
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